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Decision  C.P.M. :   Network  Reduction  and  Solution 


The  Decision  CPM  problem  may  be  termed  a  discrete  time-cost  tradeoff 
problem  generalized  to  include  the  possibility  of  interdependency  con- 
straints between  the  discrete  job  alternatives.   It  has  been  formulated 
as  an  integer  programming  problem  and  solved  by  integer  programming  and 
by  heuristic  techniques.   This  paper  shows  that  a  variant  of  the  usual 
critical  path  method  may  be  used  to  reduce  the  DCPM  networks  to  an  equiv- 
alent network  containing  only  decision  nodes  and  maximal  distances  between 
them.   This  operation,  along  with  feasibility  tests  on  the  interdependency 
constraints,  allows  substantial  elimination  of  the  precedence  constraints 
in  the  programming  formulation.   A  branch  and  bound  method  is  then  devel- 
oped to  solve  the  reduced  problem. 


Decision  C.P.M.;  Network  Reduction  and  Solution 

The  usual  network  of  the  Critical  Path  Method  has  been  generalized 

by  Crowston  and  Thompson  [2]  to  include  the  possibility  of  performing 

some  of  the  jobs  in  several  alternate  ways.   They  define  J  =  {S^ ,  S„, 

S„,..,  S  }  as  a  set  of  job  sets  comprising  a  project.   Some  job  sets  are 

unit  sets  S.  =  {S .  }   and  others  are  multi-iob  or  "decision  job"  sets 
1    ^  il  ^  J  J 

S.  =  fS.  _,  S.  „....S,  ,  /.s).   Exactly  one  job  from  each  set  must  be  per- 
1    '^  1,2'   1,2'   '  i,k(i)^         J  J 

formed  in  order  to  complete  the  project.   Associated  with  each  multi-job 

set  S.,    are  k(i)  variables  d.  ^....,    d.  ,  -..  with  the  property 
1  1,1       i,k(x) 


d. 


1  if  job  S..  is  performed 
ij    ^ 


'  0  otherwise. 

Associated  with  every  unit  set  and  multi-set  job,  S...  is  a  cost  c...  a 

ij  ij 

performance  time  t...  and  a  start  time  w...   If  it  is  assumed  that  the 
ij  ij 

project  has  a  due  date,  DD,  with  a  premium  of  "r"  dollars  per  time  period 

for  early  completion  and  "p"  dollars  per  time  period  penalty  for  late 

completion,  we  may  formulate  the  following  discrete  time-cost  tradeoff 

problem.   Find  D  =  {d..}  so  as  to  minimize  the  sum  of  job  cost  and  due 
ij 

date  penalty  or  premium,  that  is: 

N   k(i) 
Y  I 

1=1   j=l 


Minimize  Z  =  Y  T         d..  c..  -  rw^  +  pw^  (1) 

ij   ij      f    ^  f 


Subject  to: 

Interdependency  Constraints 
(a)  Mutually  exclusive; 
k(i) 

r 
j=i 


d^j  =  1       i  =  1,...,  N  (2) 
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(b)  Other;  as  required  eg. 

d.  .  <  d 
xj  -  mn 

(3) 

d.  .  =  d      etc 
ij    mn 

Precedence  Constraints 

w  +  t   <  w  (4) 

m    m  —  n 

where  S   is  an  immediate  predecessor  of  S   and  S   is  not  a  decision  job, 
m  n      m 


-  M(l  -  d_,  .)  +  w.  .  +  t.  .  <  w  (5) 

ij      iJ     ij  -  n 

where  S. .  is  a  decision  iob  and  M  is  a  large  positive  number. 

Project  Completion  Constraint 

w   -  w  +  w"  -  DD  =  0  (6) 

where  w   is  the  start  time  of  S  ,  the  last  job  in  the  project. 


In  this  formulation  a  precedence  constraint  exists  for  each  immediate 
predecessor  relationship  in  the  project  network,  and  in  total  they  con- 
strain w   to  begin  after  all  jobs  in  the  project  are  completed.   An  alter- 
nate method  of  constraining  w   is  to  allow  one  constraint  for  each  path 
from  S   to  S   in  the  network.   For  the  example  of  Figure  1,  one  of  the 
forty-five  such  constraints  would  be 

'l^   h^   ho^  hi-"   ^2.1-^^^  -  ^12, 1^  i"f 

In  this  paper  It  is  shown  that  in  the  "path"  formulation  of  the  pre- 
cedence constraints,  many  constraints  are  dominated  and  may  be  elimlnateed. 
It  is  also  shown  that  the  "other"  interdependencv  constraints  may  make  some 
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S  =  Start 


S   =  Finish 


"other"  intardependency  constraints 
S 


g^j^   only  if  S^  2  S^  ^   ^^  a^d  only  if  S 


12,2 


FIGURE  1 
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paths  infeasible  and  the  constraints  derived  from  the  paths  unnecessary. 
An  algorithm  is  developed  to  reduce  decision  networks  to  networks  that 
contain  only  decision  jobs  but  that  are  equivalent  to  the  original  net- 
work in  the  sense  that  given  any  feasible  set  of  decisions  the  resulting 
early  start  times  for  all  decision  jobs  in  the  reduced  network  will  be 
the  same  as  those  in  the  original  network.   The  reduced  network  is  then 
used  to  generate  a  set  of  undominated  path  constraints  for  the  integer 
programming  formulation.   Finally  a  branch  and  bound  algorithm  is  pre- 
sented for  the  solution  of  the  decision  network  problem. 

Dominance  Tests  for  Constraint  Elimination 

Given  a  technologically  ordered  list  of  jobs  in  the  network,  label 
the  decision  jobs 

u(l),  u(2),.,u(i),.,u(K) 

where  K  is  the  total  number  of  decision  jobs  from  all  multi-job  sets  in 

the  project.   Now  define  P  as  the  set  of  all  paths  from  S   to  S   in  the 

network.   P  may  be  partitioned  into  the  following  subsets, 
p    f  0    1    2       j       bi 

such  that 

P   =  lO  order  paths,  those  containing  none  of  the  decision  jobs) 

P   =  '1st  order  paths,  those  containing  one  of  the  decision  jobs) 

p   =  'jth  order  paths,  those  containing  exactly  j^  of  the  decision  jobs) 

p   =  {bth  order  paths,  those  containing  exactly  b^  of  the  decision  jobs} 

where  b  is  the  maximum  number  of  decision  jobs  on  any  path  from  S   to  S^, 

s     f 

It  is  now  possible  to  sub-partition  p   into  K  subsets  in  the  follow- 
ing way. 
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P^  =  IP^     P^         P^         P^    } 

where  P  ,,.    is  that  subset  of  paths  containing  only  one  decision  job 
which  is  S  ^-.v,  which  may  be  empty.   Similarily,  the  2nd  order  paths  may 
be  partitioned  into 

K! 


(K-2)!2! 


subsets  as  follows. 


^'  =  <i),u(J)>  i  =  l,2....,K-l 

j  =  i  +  1,....K 


and  P  may  be  partitioned  into 


subsets . 


(K-j)!j! 

The  maximum  number  of  subsets,  each  containing  a  unique  combination  of 
decision  jobs  will  be 

Jo  0^^ 

The  actual  number  of  subsets  required  would  be  much  less  than  this  since 
we  are  constrained  to  accept  only  one  decision  job  from  each  multi-job 
set  and  since  for  usual  projects   b  <  K. 

Now  label  each  subset  of  P-'  as  Pj^    k  =  1,2,...,     K!       and  the 

.  (K-j)!j! 

decision  jobs  contained  in  Pj^  as  Uj|  .   For  every  DCPM  problem,  any  solu- 
tion will  have  one  of  the  following  mutually  exclusive  properties.   Either 
all  j  decision  jobs  in  Uj^  will  be  performed,  or  at  least  one  of  them  is 
not  performed.   If  one  of  them  is  not  performed,  then  the  precedence  con- 
straints representing  all  paths  in  the  set  will  not  be  binding  on  w  ,  the 
start  time  of  the  finish  node.   That  is,  if 
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/  .  .   d   <  j  then  the  constraint  representing 


J    mn 


mn  for  S   €  U 
mn     K. 

any  path  in  the  subset  will  meet  the  condition 

i           2  .  -  M(l-d      )    +   t      l+T<0<w^ 

)           -        „  ,j  mn            mn I                             f 

^  mn   for   S        e  U,  ' 

mn  k 

where  T  represents  the  sum  of  the  times  of  the  unit  set  jobs  on  the  path. 
If  all  the  jobs  in  the  subset  are  performed,  then  the  finish  day  of 
the  project,  w  ,  will  be  constrained  by  equations  representing  all  paths 
in  the  subset.   Since  the  constraints  are  identical  except  that  the  value 
for  T,  the  sum  of  the  times  of  the  non-decision  jobs  on  a  particular  path, 
differ,  it  is  only  necessary  to  retain  the  constraint  that  has  the  maximum 
value  for  T.   All  constraints  representing  other  paths  in  the  subset  will 
be  dominated  by  the  one  chosen.   Thus,  we  require  only  one  path  for  each 
possible  combination  of  decision  nodes.   We  will  term  this  path  an  "un- 
dominated"  path.   It  should  be  noted  that  the  set  of  decision  jobs  Uj|  may 
later  prove  to  be  infeasible  as  the  result  of  the  "other"  interdependency 
constraints. 

Implementation  of  Dominance  Tests  for  Path  Elimination 

This  section  will  develop  an  algorithm  based  on  the  repetitive  ap- 
plication of  the  longest  path  calculation  to  generate  an  undominated  set 
of  paths  in  a  decision  network.   As  stated  above,  the  undominated  path 
will  simply  be  the  longest  path  in  a  subset  of  paths  which  all  contain  a 
particular  combination  of  decision  jobs.   To  determine  this  path  for  the 
subset  of  paths  Pj^ ,  it  is  only  necessary  to  list  the  decision  jobs,  Uj^ 
in  technological  order  and  compute  the  longest  partial  path  between  ad- 
jacent decision  jobs,  which  contains  no  other  decision  job.   The  partial  paths 
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taken  together  will  form  the  "undominated"  path  for  the  set  P^ . 

An  algorithm  based  on  the  usual  longest  path  calculation  (CPM)  was 
written  to  generate  the  maximal  distances  between  all  decision  nodes  and 

S  ,  S^.   The  algorithm  uses  two  time  values  for  each  iob  in  the  network. 

s    f  -^ 

The  "actual  time"  is  the  estimated  completion  time  for  the  job,  t   ,  which 
remains  constant  throughout  while  the  "current  time"  is  equal  either  to  actual 
time  or  -M,  a  large  negative  number.   As  we  shall  see,  if  a  job  time  is  taken 
as  -M,  then  no  path  through  the  job  has  a  non-negative  length  and  all  such 
paths  will  be  ignored.   A  flow  chart  for  the  algorithm  will  now  be  presented 
and  reference  will  be  made  to  Tables  1,  2  and  3  which  show  some  steps  of  the 
algorithm  for  the  decision  network  of  Figure  1. 

1.  List  the  jobs  of  the  decision  network  in  technological  order. 

2.  Form  a  matrix  with  a  row  for  each  job  in  the  network.   Column  1  will 
contain  the  job  number;  Column  2,  the  actual  time;  Column  3,  the 
job's  predecessors;  Columns  4,  5  and  6,  the  current  time,  the  early 
start  time  and  the  early  finish  time  (i.e.  Table  1). 

3.  Identify  all  decision  jobs  in  Column  1.   (in  Table  1  brackets  are 
used) . 

4.  Set  the  current  time,  Column  4,  for  all  decision  jobs  at  -M,  for  all 
other  jobs  set  the  current  time  equal  to  actual  time.   Set  counter 

n  =  S  .   Go  to  Step  8. 

5.  Beginning  at  H,  search  down  Column  1  for  the  first  decision  job  be- 
low n,  say  this  is  S.  ..   Set  11  =  S.  .. 

6.  Set  the  early  start.  Column  5,  for  all  jobs  above  n  at  -M, 

7.  Set  the  current  time  of  H  and  the  early  start  of  n  at  0. 

8.  Calculate  the  early  start  time  of  all  jobs  beneath  n  on  the  tech- 
nologically ordered  job  list.   If  the  early  start  value  calculated 
is  negative,  enter  -M  in  the  early  start  column.   To  calculate  the 
early  start  of  the  job  directly  beneath  11,  list  all  the  predecessors 
of  the  job  and  for  each  predecessor  add  the  predecessor's  early  time 
to  its  current  time.   The  largest  of  these  values,  for  all  predeces- 
sors, is  the  early  start  time  for  the  job  being  considered.   Then 
proceed  to  the  next  job  on  the  technologically  ordered  list  for  the 
next  early  start  calculation, 

9.  For  all  decision  jobs,  add  early  start  time  to  actual  time  to  com- 
pute early  finish  time.  Column  6. 

10.   For  all  decision  jobs  beneath  n  on  the  technological  list  and  for 

S  ,  record  the  early  finish  time  if  it  is  non-negative.   These  times 
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will  be  the  longest  sub-path  from  the  finish  of  decision  job  n  to  the 
finish  of  the  decision  job  being  examined  (the  longest  path  contain- 
ing no  other  decision  job) . 

11.  If  n  is  the  last  decision  job  in  the  technological  list,  go  to  Step 
12.   If  not,  go  to  Step  5. 

12.  HALT. 

This  algorithm  has  been  applied  to  the  problem  of  Figure  1  with  results 
shown  in  Table  3.   Each  non-blank  entry  in  this  matrix  indicates  the  length 
of  the  maximal  path  between  pairs  of  decision  jobs  or  between  the  Start  and 
Finish  nodes.   Note  that  the  entries  for  line  one,  the  maximal  distances 
from  S   to  all  decision  jobs,  are  taken  directly  from  Column  6,  Table  1. 
Similarly,  the  entries  for  row  two,  the  distances  from  S,  ,  to  all  decision 

0.1 

jobs,  are  taken  from  Column  6  of  Table  2.   Blank  entries  in  the  matrix  indi- 
cate that  no  path  not  containing  a  decision  job  connects  the  two  jobs. 
Figure  2  shows  the  "reduced"  decision  network  derived  directly  from  the  matrix 
of  Table  3. 


Feasibility  Tests  for  Path  Elimination 

Consider  a  particular  subset  of  paths  P^.   For  the  constraints  rep- 
resenting these  paths  to  constrain  w  ,  the  following  condition  must  hold. 

S  d   =  j 

mn  for  S   £  U, 
mn    k 

This  condition  may  be  contradicted  by  the  "other"  interdependency  con- 
straints, so  that  the  set  U.-'  is  infeasible  and  the  undominated  path  from 

k 

Pj^  may  be  eliminated  from  the  integer  programming  formulation  of  the  DCPM 
problem.  Paths  not  eliminated  by  the  interdependency  constraints  will  be 
termed  "feasible"  paths. 
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Application  of  Dominance  and  Feasibility  Tests  to  a  DCPM  Network 

When  the  tests  of  the  previous  three  sections  are  applied  to  precedence 
equations  of  the  example  Figure  1,  the  reduced  set  of  equations  of  Table  4 
result.   For  example,  the  set  P/-.^  ,x  contains  the  6  paths 

•^1^  ^7  -'ho  -"(^-'^12,1^  -^^2.1^  ^6  I'^f 

h  ^  "^7  ^  ho  ^  hi   -^(l-'^U.P  ^  ^2,1  ^  ^6  ^w^ 

h'-h^   ^8^  ho   +  '^ll  -^(1-^12.1^  ^  ^2,1  +  t,,  <  w, 

lb  —  1 

h^  h^  h^  ho  -^<i-^i2.i>  -^  ^2.1  -^  he  l"f 
h^h^  '8^  ho  ^  hi  -"(^-^12.1^  "■  'i2.i  ■"  he  l«f 

h  -^  '^A  ^  ^8  -^  ho  -^(^-^12. 1>  -^  h2.1  ^  hb  l^f 
These  may  be  rewritten 

-M(l-d^2  i>   +  38  £  w^ 

-M(l-d^2  i>  +  ^3  £w^ 

-M(l-d^2  i>  +  ^6  iw^ 

-M(l-d^2  !>  +  ^1  1  w^ 

-M(l-d^2  i)  +  ^A  £  w^ 

-M(l-d^2  1^  +  39  £  w^ 

With  -M(l-d^„  ^)  +  46  <_ w  ±n  the  problem,  the  other  5  constraints  are  re- 
dundant. 


Calculation  of  Maximal  Distances  from  S   to  all  Decision  Jobs 

s 

12                 3  4        5       6 

Job     Actual         Predecessors  Current   Early   Early 

Time  Time     Start    Finish 

0  0 

^s  12 

S  10 
s 

2,3  4 

3  4 

1  -M 
1  -M 
1  -M 

1,(6,1)  2 

4,(6,3)  3 

5,8  -M 

5  -M 

5  -M 

(6, 2), (6, 3), 7, 8  10 

10  5 

10,11  -M 

10  -M 

8,(9,1)  2 

(9, 2), (9, 3), 13  8 

Sg,(9,3)  -M 

Sg  ,  (9  ,  3)  -M 

(12,1), (12, 2), 14, (15,1)  10 

Sg,(15,2)  -M 

Sg,(15,2)  -M 

S^, (17,1), (17, 2)  5 

11,16,18  0 

Table  1 

(Decision  Job) 


s 

s 

1 

0 
12 

2 

10 

3 

8 

4 

4 

5 

4 

(6,1) 

4 

(6,2) 

6 

(6,3) 

8 

7 

2 

8 

3 

(9,1) 

5 

(9,2) 

10 

(9,3) 

15 

10 

10 

11 

5 

(12,1) 

4 

(12,2) 

5 

13 

2 

14 

8 

(15,1) 

9 

(15,2) 

13 

16 

10 

(17,1) 

11 

(17,2) 

3 

18 

5 

S, 

0 

0 

0 

0 

12 

0 

10 

0 

8 

10 

14 

8 

12 

12 

16 

12 

18 

12 

20 

12 

14 

14 

12 

17 

22 

12 

22 

12 

27 

17 

27 

27 

32 

32 

36 

27 

32 

17 

19 

19 

27 

0 

9 

0 

13 

27 

37 

0 

11 

0 

3 

0 

5 

37 

37 
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Calculation  of  Maximal  Distances  From  S,  ,  to  All 

6,1 

Successor  Decision  Jobs 
12  3  4       5        6 

Job    Actual  Predecessors         Current  Early    Early 

Time  Time   Start    Finish 

0 

S 
s 

S 
s 

S 
s 

2,3 
3 
1 

1  -M  -M  -M 

1  -M  -M  -M 
1,(6,1) 

4,(6,3)  3  -M  -M 

5,8  -M  -M  -M 

5  -M  -M  -M 

5  -M  -M  -M 

(6, 2), (6, 3), 7, 8  10  2  12 

10  5  12  17 

10,11  -M  17  21 

10  -M  12  17 

8,(9,1)  2  -M  -M 

(9, 2), (9, 3), 13  8  -M  -M 

S  ,(9,3)  -M  -M  -M 

S  ,(9,3)  -M  -M  -M 

(12,1), (12, 2), 14, (15,1)  10  -M  -M 

S  ,(15,2)  -M  -M  -M 

S  ,(15,2)  -M  -M  -M 

5^,(17,1), (17. 2)  5  -M  -M 

11,16,18  0  17  17 


Table  2 


< 

5 

s 

0 

1 

12 

2 

10 

3 

8 

4 

4 

5 

4 

(6, 

,1) 

4 

(6: 

,2) 

6 

(6, 

.3) 

8 

7 

2 

8 

3 

(9: 

,1) 

15 

(9, 

,2) 

10 

(9: 

.3) 

15 

LO 

10 

LI 

5 

(12: 

,1) 

4 

(12, 

,2) 

5 

13 

2 

14 

8 

(15, 

,1) 

9 

(15, 

,2) 

13 

16 

10 

(17, 

,1) 

11 

(17, 

,2) 

3 

18 

5 

S 

0 

0 

-M 

12 

-M 

10 

-M 

88 

-M. 

4 

-M 

4 

-M 

0 

0 

-M 

-M 

-M 

-M 

2 

0 

3 

-M 

-M 

-M 

-M 

-M 

-M 

-M 

10 

2 

5 

12 

-M 

17 

-M 

12 

2 

-M 

8 

-M 

-M 

-M 

-M 

-M 

10 

-M 

-M 

-M 

-M 

-M 

5 

-M 

0 

17 
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Reduced  Network  Matrix 


T  0 


! 
i 

6,1 

6,2 

6,3 

9,1 

9,2 

9,3 

12,1 

12,2 

15,1 

15,2 

17,1 

17,2 

h 

'  s 

i   s 

16 

18 

20 

22 

22 

27 

36 

32 

9 

13 

11 

3 

37 

■  ^'^ 

21 

17 

17 

;  6,2 

19 

15 

15 

6,3 

8 

22 

18 

23 

;  9,1 

20 

'  9,2 

18 

9,3 

9 

13 

18 

'12,1 

10 

12,2 

10 

!15,1 

10 

j  15,2 

11 

3 

17,1 

5 

1  17,2 

5 

Table   3 
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FEASIBLE  NETWORK  PATHS 


-M(l-d,  ,)  +  33  <  w^ 

D,i  —   f 

-M(l-d,  „)  +  43  <  w^ 

o ,  J  —  I 


-M(l-dg  ^)  +42  £w^ 
-M(l-dg  2)  +  40  £w^ 
-M(l-dg  2)  +  45  £w^ 

-M(l-d^2  1^  +  46  £w^ 
-M(l-d^^  ^)  +  19  ^w^ 
-M(l-d^^  ^)  +  16  £w^ 
-M(l-d^^  2^  "^  ^  -"^f 
-M(l-d^  ^)  -M(l-d^2  1^  +  47  £w^ 
-M(l-d^  2>  -M(l-d3^2  2^  +  43  <  w^ 
-M(l-d^  ^)  -M(l-d^2  1^  +  49  ^w^ 
-M(l-dg  ^)  -M(l-d^^  j^)  +  46  £w^ 

-M(l-d^3^2>  -^(1-^17, 1>  ■^29^w^ 
-M(l-d^5  2^  -M(l-d^^  2>  +  21  £w^ 
-Md-dg^^)  -M(l-d^5^2^  -M(l-d^^  ^)  +  56  £  w^ 
-M(l-dg^3)  -M(l-d^^^2)  -M(l-dj^7  2^  +  48  £  w^ 


Table  4 
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We  will  now  give  an  example  of  constraint  elimination  by  feasibility 

tests. 

The  inequality 

-M(l-d,  -)  -M(l-d    )  +  43  <  w  will  be  potentially  binding 
6  J 1         -L^  >  ^  t 

only  if 

'^e.i'' ^^12,2  =  2 

However,  the  interdependency  constraints 

^6,2  =  ^12.2    ^""^       ^6,1  +  ^6,2  +  ^6.3  =  ^ 
together  violate  this  condition.   Therefore,  the  path  constraint  may  be 
dropped. 

Lower  Bound  Calculation 

The  set  of  path  constraints  has  now  been  reduced  from  one  per  path  in 

the  original  decision  graph  to  one  per  combination  of  decision  jobs  lying 

in  a  path  from  S   to  S^.   The  number  may  be  even  further  reduced  because  of 
s     f 

feasibility  tests  on  the  interdependency  relations. 

We  will  now  define  set 
M  =  {all  feasible,  undominated  paths} 
and  set 

M.  =  {the  subset  of  paths  in  M  that  contain  any  of  the  decision  jobs  in  de- 
cision set  S  ,  a  multi-job  set} 
Each  path  in  M.  contains  exactly  one  of  S   ,    j  =  1,2,...  k(i)  since  we  are 
assuming  that  only  one  job  from  each  multi-job  set  may  be  selected,  but  they 
may  also  contain  other  decision  jobs. 

We  will  now  define  a  "feasible  solution"  to  a  DCPM  problem  as  a  selec- 
tion of  a  set  of  decision  jobs  that  will  satisfy  all  interdependency 
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constraints  on  the  problem.   This  set  of  decision  jobs  will,  with  the  non- 
decision  jobs,  form  a  usual  project  graph  of  the  critical  path  method.   Given 
the  project  graph  resulting  from  a  "feasible  solution,"  we  state  that  one  of 
the  paths  from  M.  will  be  included  in  that  graph.   If  no  such  path  exists, 
then  no  path  through  S.  exists  in  the  network.   This  would  violate  the  assump- 
tion that  S   is  a  predecessor  and  S   is  a  successor  to  all  jobs  in  the  project 
graph. 

Each  constraint  representing  a  path  from  M.  will  be  of  the  general  form 

-M(l-d   ) -M(l-d   )  +  T    <  w 

ij  mn     i,e  —  t 

where  T.    is  the  length  of  path  e  EM.. 

If  one  of  the  paths  from  M.  must  be  performed,  then  a  lower  bound  on 

w-  is  the  lowest  value  of  T.    for  all  paths  in  M. .   Furthermore,  if  we  cal- 
f  i,e         ^         1 

culate  a  lower  bound  for  every  decision  node,  we  may  select  the  highest  of 
the  lower  bounds  as  a  lower  bound  on  w  .   Any  undominated,  feasible  path  which 
is  shorter  than  this  bound  may  be  eliminated  since  it  is  redundant,  given  the 
bound  on  w^. 


An  Application  of  Lower  Bound  Calculation 

The  reduced  set  of  equations  from  Table  4  will  be  used  to  illustrate 
the  dominance  tests  described  above.   The  only  feasible  paths  through  S,, 


that  is  the  set  M  ,  are  the  following  paths, 


-n(l-d,  ,)  +  33  <  w^ 
0  ,  i       —  t 

-M(l-d,  „)  +  43  <  w^ 
D  ,  J       —  r 

-M(l-d^  ^)  -M(l-d^2  !>  +  ^7  £w^ 
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COMPUTATION  OF  LOWER  BOUND  ON  w. 


p 

6,1  6,2 

6,3 

9,1  9,2  9,3 

12,1  12,2 

15,1  15,2 

17,1  17,2 

1 

f 

o 

p 

37 

1 
p 

1 

1 

1 

1 

1 

1 

1 

1 

1 

33 

43 
42 

!  40 

45 

i 

46 

19 

! 

16 
8 

2 
P 

1 

1 

1 

1 

1 

1 
1 

1 

1 
1 

1 

1 

47 
43 
.9 

1 

! 

29 

1 
21 

3 
P 

1 

1 

1 

'  56 

1 

1 

1 

1  48 

BOUND 

33 

40 

43 

19 

8 

Table  5 
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-M(l-d^  ^)    -M(l-d^2  2^  +  ^3  ^w^ 
-M(l-d^  2>  -M(l-d^2  i>  +  A9  ^w^ 


The  minimum  of  the  T.      e  =  1, . . .  ,4  is  33  and  therefore  33  is  the 
i,e 

lower  bound  provided  by  M, .   The  bound  provided  by  M   is  40,  M   is  43, 
M   is  19,  and  M   is  8.   These  calculations  are  shown  in  Table  5.   The 
highest  of  the  bounds  is  then  43,  and  all  paths  shorter  than  or  equal  to 
this  may  be  removed.   For  this  particular  problem  this  reduction  technique 
does  not  give  any  further  improvement.   The  final  set  of  predecessor  con- 
straints are  shown  in  Table  6.   The  equations  (1)  -  (6)  with  the  reduced 
set  of  precedence  constraints  may  now  be  solved  by  usual  integer  program- 
ming methods. 


REDUCED  SET  OF  PRECEDENCE  CONSTRAINTS 


+  43  >^  w 

-M(l-dg  ^)  +  45  ^  w^ 

-M(l-d^2  i>  +  46  ^  w^ 

-M(l-dg^3)  -M(l-d^5^,)  +46  >.w^ 

-M(l-d^^^)  -M(l-d^2,l>  +^7^w^ 

-M(l-dg  3)  -M(l-d^5  2)  -Md-dj^y  2^  +  48  ^w^ 

-M(l-d^  3)  -M(l-d^2  !>  +  ^9  >^  w^ 

-M(l-dg  3)  -M(l-d^^  2>  -M^l-dj^y  {)  +  56  ^  w^ 


Table  6 
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Branch  and  Bound  Algorithm 

We  now  present  a  branch  and  bound  algorithm  for  the  problem  stated 
in  equations  (1)  -  (6).   The  word  "branch"  in  this  context  refers  to  a 
specific  decsiion  rule  for  implicitly  enumerating  a  tree  of  all  feasible 
solutions  to  a  problem.   The  "bound"  term  refers  to  the  computation  of 
a  lower  bound  (minimization  problem)  on  the  value  of  the  criterion  func- 
tion at  each  node  of  the  enumeration  tree.   This  algorithm  will  be  based 
on  the  reduced  set  of  constraints  of  the  previous  section. 

The  bounding  calculation  is  based  on  the  decomposition  of  the  full 
problem  into  two  related  sub-problems  and  on  the  factoring  of  the  K  de- 
cision jobs  into  three  mutually  exclusive  subsets, 

SE  =   the  subset  of  decision  jobs  that  have  been  excluded  from  any 
solution  D 

SA  =   the  subset  of  decision  jobs  that  have  been  accepted  for  in- 
clusion in  D 

SR  =  all  remaining  decision  jobs 

The    first   sub-problem  solves    for   a   lower  bound   C.    on    the   cost   of   the 

J 

decision  jobs  at  a  node  in  the  enumeration  tree,  given  all  interdepen- 
dency  constraints  and  SE.   That  is 

Minimize 

N   k(i) 
C   =   Z    i:     c.  .  d.  .  (7) 

'        i=l  j=l    ^J   ^J 

Subject  to 

Interdependency  Constraints 

Mutually    kCi) 

Exclusive     I   d..  =  1    i  =  1,...,N  (8) 

j  =  l   '' 

Other       d. ^  =  d    etc.  (9) 

ij    mn 

Exclusion  Constraints 

d. .  =  0         d. .  e  SE  (10) 
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If  no  "other"  interdependency  constraints  exist,  the  problem  is 
solved  trivially  by  choosing  for  each  multi-set  job  S..  ^  SE.   For 
problems  with  certain  "other"  constraints  the  modified  shortest  path 
algorithm  of  Crowston  [1]  will  efficiently  compute  C  .   In  general,  the 
problem  of  equations  (7)  -  (10)  may  be  solved  as  a  linear  program  for  a 
lower  bound  on  C   or  a  0-1  tree  search  algorithm  could  be  developed. 

The  second  sub-problem  is  to  solve  for  a  lower  bound,  C  ,  on  the 
cost  of  the  due  date  penalty  (or  reward). 

Minimize  C   =  -rw  +  pw  (11) 

Subject  to: 

Precedence  Constraints 


{I 


mn  for  S   c   U? 
mn    k 


-M(l-d   )  +  t   >  +  T  <  w^  (12) 

mn     mn  ;     —  f 


for  all  feasible  paths 

Acceptance  Constraints 


d.  .  =  1  S. .  e  SA  (13) 


Subsequent  discussion  of  the  branching  strategy  will  show  the 

development  of  the  set  of  decision  jobs  SE.   Given  SE  and  the  mutually 

k(i) 

exclusive  interdependency  constraints  of  the  multi-iob  sets,  S      d,.  =  1, 

j  =  l   ^1 

then  the  inclusion  of  some  jobs  in  SA  becomes  mandatory.   Given  SA  and 

either  the  reduced  network  of  Figure  2  or  the  reduced  equation  set  of 

Table  6,  a  lower  bound  on  W,  and  thus  a  bound  on  C   can  be  simply  deter- 

f  p  ^  ■' 

mined.   For  example,  if  SE  =  | S^  ^,  S^  ^,    S^^  2  I  '  ^hen  SA  =  {  S^  ^,  S^^  ^f 
and  from  Table  6  a  lower  bound  on  W   is  46  days  and 
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.   p(46-DD),  46  >  DD 
C   =  j 
P    (  -r(DD-46)  otherwise 

Thus  the  lower  bound  on  total  cost,  C  ,  for  the  partial  solution  is  Cj  + 

C  .   In  addition,  if  it  is  found  that  the  complete  solution  D  to  the  job 
P 

cost  minimization  problem  also  gives  a  W  resulting  in  cost  C  ,  we  may 
terminate  further  search  down  the  path.   In  this  case,  the  optimum  solution 
containing  the  partial  solution  has  been  found  since  we  have  a  complete  so- 
lution with  minimum  job  cost  and  path  length,  given  SE,  the  general  branch- 
ing strategy  to  be  used  here  is  as  follows.   From  the  set  of  all  enumeration 
paths,  the  set  \jith  the  maximum  number  of  excluded  jobs  will  be  taken.   For 
each  of  these  paths  a  minimum  bound  will  be  calculated  and  the  path  with 
the  minimum  lower  bound  will  be  selected  for  further  elaboration,  unless 
it  can  be  shown  that  a  "complete"  solution  exists.   Once  the  path  has  been 
chosen,  a  decision  rule  is  required  to  select  decision  jobs  that  will  be 
examined  at  the  next  stage.   The  approach  to  be  used  here  is  similar  to 
that  used  by  Eastman[4]  and  later  by  Shapiro  [5]  for  the  travelling  sales- 
man problem.   The  computation  of  C   at  a  node  implies  the  selection  of  a 
feasible  set  of  decision  jobs  at  that  point  in  the  tree.   Given  this  set 
of  jobs  it  is  again  possible  to  use  the  reduced  set  of  precedence  con- 
straints to  find  the  length  of  the  implied  critical  path  and  also  to  iden- 
tify the  decision  jobs  on  the  critical  path.   Note  that  the  critical  path 
determined  here  will  ordinarily  not  be  the  same  as  that  found  in  the  cal- 
culation of  C  ,  which  was  based  solely  on  those  jobs  S, .  £  SA. 
p'  ^  -J      ij 

The  critical  path  variables  associated  with  C  then  form  one  branch 
from  the  new  node,  the  other  branches  being  formed  by  each  variable,  in- 
dividually excluded.   Any  solution  to  the  problem  must  contain  either  the 
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set  of  all  decision  nodes  on  the  original  critical  path,  or  have  at  least 
one  of  then  excluded.  If,  for  example,  at  a  particular  node  on  the  tree, 
the  cheapest  feasible  set  of  decisions  contained  S    and  S     and  both 

D  ,  i  1^  ,  1 

of  these  jobs  were  on  the  critical  path,  we  would  branch  as  shown  below. 


In  the  application  of  the  idea  to  the  DCPM  problem,  the  combination  of 

jobs  on  the  critical  path,  say  S,  ,,  S^^  ,  is  feasible  rather  than  infeas- 

0,1   iZ ,  1 

ible  as  in  Shapiro's  problem.   However,  although  it  is  feasible,  it  is 
never  necessary  to  develop  that  particular  node  further.   Given  the  ex- 
clusions and  acceptances  to  that  point  on  the  path,  a  programming  problem 

is  solved  to  minimize  the  total  iob  cost,  C,.   The  minimum  cost  feasible 

-"  J 

set  of  jobs  was  then  tested  in  the  equation  of  Table  6  to  determine  pro- 
ject length  and  project  length  cost  C  .   At  this  point,  the  bound  C   = 

C,  +  C   is  also  the  true  cost  of  a  feasible  solution  to  the  problem.   Thus 
J    p 

the  node  is  completely  developed. 

If  at  any  stage  of  the  development  of  the  tree,  ail  new  paths  gener- 
ated are  bounded  by  a  previous  complete  solution  or  are  infeasible  or  if 
a  complete  solution  has  been  obtained,  it  is  necessary  to  move  backwards 
up  the  tree  and  pick  new  paths  to  develop.   Again  we  find  the  unbounded 
paths  that  have  the  maximum  number  of  job  acceptances  and  exclusions  and 
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of  these  select  the  one  with  the  minimum  bound.   This  implies  that  we  move 

up  the  tree  to  the  first  open  node.   The  algorithm  will  now  be  presented 

in  flow  chart  form,  then  applied  to  the  sample  problem  of  Figure  1. 

Step  1.   Reduce  the  network  by  the  dominance,  feasibility  and  lower 
bound  tests  to  obtain  the  reduced  constraint  matrix. 

Step  2a.   Solve  the  problem  of  equations  (7)  -  (10)  to  find  the  min- 
imum cost  set  of  decision  jobs  that  will  meet  all  interdependency 

constraints.   Set  the  value  of  C   to  equal  the  cost  of  this  solution. 

P 
Go  to  Step  3. 

Step  2b.   Select  the  complete  solution  and  the  job  cost  C  calculated 

at  Step  7  for  the  path  to  be  elaborated. 

Step  3.   Calculate  the  length  of  the  critical  path,  w  ,  associated 

with  the  decision  jobs  selected  at  Step  2.   This  may  be  determined 

either  from  the  original  project  graph  or  from  the  reduced  set  of 

constraints.  Table  6.   Let  C  =  -rw   +  pw^  given  w  -  w^  "*"  ^f  =  DD. 

Thus,  C   is  wither  the  early  premium  or  late  penalty  associated  with 

finish  day  w  and  due  date  DD.   Now  establish  one  branch  from  the 

first  node  which  includes  the  complete  solution  of  Step  2,  that  is 

a  particular  set  of  decision  jobs,  with  a  total  project  cost  of 

C  +  C-  =  C  . 
T    J     p 

Step  4.   Establish  alternate  branches  from  the  current  node  by  adding 
nodes  which  specifically  exclude  each  of  the  decision  jobs  on  the 
critical  path  determined  at  Step  3. 

Step  5.   For  each  new  branch  update  SE,  then  solve  the  programming 
problems  of  equations  (7)  -  (10)  for  C  . 

Step  6.   For  each  new  branch  use  the  updated  set  SE  and  the  mutually 
exclusive  interdependency  constraints  to  update  SA.   Solve  the  prob- 
lem of  equations  (11)  -  (13)  for  a  lower  bound  on  C   using  the  re- 
duced set  of  precedence  constraints  (Table  6).   Calculate  a  bound  on 

total  cost  C_  =  C ^  +  C   for  each  new  node. 
1    J     p 

Step  7.   Record  the  full  solution  to  the  programming  problem  of  Step 
5  and  calculate  the  critical  path  length  associated  with  each  solu- 
tion. 
Step  8.   Test  all  partial  solutions  from  the  last  node  for  feasibility^ 
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complete  solution  or  bounded  solution. 

(a)  If  the  program  of  Step  5  has  no  solution,  then  the  partial 
solution  being  tested  is  infeasible  and  no  further  search  on 
that  partial  solution  is  required. 

(b)  If  for  any  path  the  bound  of  Step  6  is  equal  to  the  complete 
solution  cost  of  Step  7,  then  we  have  a  complete  solution  to 
the  problem  and  no  lower  cost  will  be  found  containing  this 
particular  partial  solution.   If  the  total  cost  of  the  full 
solution  is  less  than  the  cost  of  the  existing  best  solution, 
update  the  existing  best  with  the  new  value. 

(c)  If  the  lower  bound  calculated  in  Step  6  is  higher  than  an 
existing  complete  solution  to  the  problem,  the  partial  solu- 
tion is  bounded  and  need  not  be  considered  further. 

Step  9.   Choose  from  the  feasible,  unbounded  partial  solutions  at  the 
current  node  the  partial  solution  with  the  minimum  cost  bound.   If  a 
tie  exists,  choose  the  path  with  the  maximum  number  of  forced  accep- 
tances.  Go  to  Step  2b. 

If  no  feasible,  unbounded  partial  solution  exists,  then  the 
node  is  closed  and  we  backtrack  one  level  in  the  graph.   Go  to  Step 
8. 

If  all  nodes  in  the  network  are  closed,  then  the  problem  is 
solved  and  the  current  best  solution  is  the  optimal  solution. 


Application  of  the  Reduced  Constraint  Algorithm 

The  steps  of  the  algorithm  will  now  be  illustrated  for  the  decision 
network  of  Figure  1  with  DD  =  45  days,  r  =  $20  and  p  =  $40.  The  complete 
solution  tree  is  shown  in  Figure  3. 

Step  1:   Construct  the  reduced  set  of  equations  of  Table  6. 
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Step  2 •   Solve  the  problem  of  equations  (7)  -  (10). 

for  Cj  =  $50.  The  solution  consists  of  S^  2'  ^g    y    ^^^  i>  ^^^  2'  ^17  i' 

Step  3:   Using  the  data  of  Table  6  we  find  that  S    ,  S     ,  and  S 

would  lie  on  a  critical  path  56  days  in  length.   C  =  (56-45)  (^iiO)  =  $440. 

C_,  =  C  +  c  =  50  +  440  =  $490. 
1    J    p 


Q 


C.P.  JOBS   OTHER  JOBS 


$490 


V^9,3'^15,2'  ^17,1 


^6,3'^12,1 


COMPLETE  SOLUTION 

Step  4:   We  now  establish  alternate  branches  excluding  all  jobs  on  the 
critical  path  of  Step  3 . 


EXCLUDE  S . . 
ij 


$490 


COMPLETE  SOLUTION 


Ste£_^:   For  each  excluded  job,  S^  ,  we  determine  the  remaining  minimum 


cost,  C   feasible  solution. 
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© 


Step  6:   For  each  exclusion  calculate  the  minimum  bound  for  the  length 
of  the  critical  path.   If  we  exclude  S^  „,  then  we  force  the  acceptance 
of  no  job,  and  therefore  the  minimum  length  of  the  critical  path  would 
be  the  absolute  minimum,  43  d^ys.   This  gives  a  reward  of  (45  -  43) (20) 
=  -$40.   If  we  reject  S    ^   „,  then  we  are  forced  to  accept  S^ -  ,  .   However, 
no  path  that  includes  only  decision  job  S^ _  ^  is  longer  than  43  days. 
Therefore,  the  minimum  completion  cost  is  again  -$40. 


P 


Step  7  ;   Compute  the  complete  cost  for  each  node  developed  in  Step  4, 
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S.  . 


9,3 


15,2 
(^15,1) 


17,1 
(^7,2> 


Solution  Given  C 


^6,3'   ^9,2'   ^12,1'   ^15,2'   ^17,1 


^6,3'   ^9,3'   ^12,1'   ^15,1'   ^17,1 


^6,3'   ^9,3'   ^12,1'   ^15,2'   ^17,2 


c 
J 

"f 

C 
P 

150 

49 

160 

150 

49 

160 

200 

49 

160 

310 


310 


360 


Step  8«   A.11  three  paths  S-  _,  S^^.  „,  S     are  feasible  and  all  but 

S  ^    are  unbounded.   The  minimum  cost  solution  available  costs  $310, 

Step  9:   Select  node  S   ^  for  elaboration  since  both  it  and  S^_  „ 

have  a  minimum  bound  of  C   =  $110. 

The  final  solution  as  shown  in  Figure  3  contains  iobs  S,  „,  S_  „, 

^  -■     6,2   9,2 

;.„  „,  S,^  „,  and  S,^  ,  with  costs,  C.  +  C  =  C^,  300  -  40  =  $260. 
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Summary  and  Conclusions 

This  paper  demonstrates  techniques  for  the  elimination  of  redun- 
dant precedence  constraints  in  the  integer  programming  formulation  of  the 
Decision  CPM  problem.   It  is  shown  that  only  one  constraint  must  be  in- 
cluded for  every  combination  of  decisions  jobs  lying  on  a  path  in  the 
project  graph.   This  number  is  further  reduced  by  testing  the  feasibil- 
ity of  combinations  of  jobs  in  each  path,  given  the  interdependency  con- 
straints on  the  problem  and  by  lower  computation  of  a  lower  bound  on  the 
project  length.   These  techniques  should  be  useful  for  problem  simplifi- 
cation even  if  heuristic  methods  are  to  be  used  in  solving  the  problem. 
Finally  a  branch  and  bound  method  is  described  which  makes  use  of  the 
reduced  set  of  precedence  constraints  to  efficiently  determine  an  optimal 
solution  to  the  problem. 
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